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ABSTRACT 



We analytically compute the orbital effects induced on the motion of a spin- 
ning particle geodesically traveling around a central rotating body by the general 
relativistic two-body spin-spin and spin-orbit leading-order interactions. Con- 
cerning the spin-orbit term, we compute the long-term variations due to the 
particle's spin by finding secular precessions for the inclination I of the particle's 
orbit, its longitude of the ascending node Q and the longitude of pericenter w. 
Moreover, we generalize the well-known Lense-Thirring precessions to a generic 
orientation of the source's angular momentum by obtaining an entirely new effect 
represented by a secular precession of J, and additional secular precessions of Q 
and vj as well. The spin-spin interaction is responsible of gravitational effects 
a la Stern-Gerlach consisting of secular precessions of /, i?, w and the mean 
anomaly M.. Such results are obtained without resorting to any approximations 
either in the particle's eccentricity e or in its inclination J; moreover, no preferred 
orientations of both the system's angular momenta are adopted. Their generality 
allows them to be applied to a variety of astronomical and astrophysical scenarios 
like, e.g., the Sun and its planets and the double pulsar PSR J0737-3039A/B. It 
turns out that the orbital precessions caused by the spin-spin and the spin-orbit 
perturbations due to the less massive body are below the current measurability 
level, especially for the solar system and the Stern-Gerlach effects. Concerning 
the solar Lense-Thirring precessions, the slight misalignment of the solar equa- 
tor with respect to the ecliptic reduces the gravitomagnetic node precession of 
Mercury down to a 0.08 milliarcseconds per century level with respect to the 
standard value of 1 milliarcseconds per century obtained by aligning the z axis 
with the Sun's angular momentum. The new inclination precession is as large as 
0.06 milliarcseconds per century, while the perihelion's rate remains substantially 
unchanged, amounting to —2 milliarcseconds per century. Further studies may 
be devoted to the extrasolar planets which exhibit a rich variety of orbital and 
rotational configurations. 

Subject headings: Classical general relativity; Experimental tests of gravitational 
theories; Celestial mechanics; Mercury; Pulsars 
PACS: 04.20.-q; 04.80.Cc; 95.10.Ce; 96.30.Dz; 97.60.Gb 
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1. Introduction 



In this paper, we deal witli tlie issue of tlie dynamics of massive spinning bodies in 
a gravitational field in the framework of the weak-field and slow-motion leading-order 
approximation of th^ General Theory of Relativity (GTR). In particular, we will 
analytically work out the long-term orbital effects experienced by an object of mass m 
endowed with proper angular momentum S in motion around a central body of rnass M 
and proper angular momentum J. We will also generalize the iLense and Thirringi (119181 ) 
effect, caused by J on the orbital motion of a moving particle, assumed non-spinning, to 
arbitrary spatial orientations of it; it may have important consequences from an empirical 
point of view when data from specific astronomical scenarios are analyzed. The case of 
arbitrary masses Ma and Mb will be taken into account as well. 

For a general review on the gr avitational spin-dependent tw' o -body pro blem, see 
Barker Sz O'Connelll ( 119791 ) ; see also iKhriplovich and PomeranskyI ( 119991 ) and iKhriplovich 
(l2008h . 



According to [Barker fc O' Conned (Il970[ ) , the leading-order spin-spin and spin-orbit 
perturbations to the reduced Hamiltonian are, in the slow-motion limit and for the extreme 
mass ratio case M ^ m. 



"Ho- 



-^[3(J.r) {cT.r)-rHcT.J)], 



3GM 



(T-L) 



In it G is the Newtonian constant of gravitation and c is the speed of light in vacuum. 
Moreover, cr = S/m and L = r X v are the reduced spin and orbital angular momentum 
of the particle moving with speed v at distance r from the central body: [a] = [L] = 
T~^, so that [T-La] = [Hah] = T~^, i.e. they have the dimensions of an energy per unit 
mass. We point out that Eq. ([1]), and also Eq. (|T3l) treated later on to generalize the 
Lense-Thirring ef fect, are a direct conseque nce of the full spin-orbit Hamiltonian for two 
arbitrary masses (IBarker fc 0'Connedll975l ): see Eq. (132|) and Eq. (IM|) -Eq. (|35|1 in Section 
13.21 for the connections between such scenarios. 



The spin-spin term "H^ wBS_obtained for the fi rst time bv ISchifj (119601). w hile the 
spin-orbit teru^'HaL is due to iFokkerl ( 19211 ): see also [Barker &: O'Connelll ( 1970l ) for another 
derivation followi ng a different line of reasoning. J-Lrr yields a gravi t ational acceleration d 



la Stern-Gerlach ( Barker fc O'Connelll 1 19 70l . Il975l : iMashhood l2000l : IChicone et all 12005 



^For a post-Newtonian, lea ding-order treatment including the PPN parameters (3 and 7, 



see 



Barker fc O'Connelll (Il976t ). 



^It is the gravitational analogue of the Thomas precession in electromagnetism which, 
curiously, was obtained 5 years later (jThomaslll926l ). 
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Faruqud |2006| ) whose eff ects on the orbital motion of a spinning particle were recently 
worked out in detail by llorid ( 120031 ) under certain simplifying assumptions about the 
orientation of J and the overall configuration of the system. For other studies concerning 
vari ous consequences of t h e gravitational St er n-Gerlac h effec ts on the orb i tal motion, see, 



e^ 



Barker fc O'Connelll fll970l. Il975l. Il979h: iFaruaud (1200 



Chicone et al.l ( l2005l ): iMashhoon and Singhl (120061): iFaruaue 



spin-orbit term were derived by lBarker fc O'Connelll (11970 



) : iBini et al.l (12004 120051 ): 
(l2006h . Orbital effects of the 
1975h . Concerning the role of 



gen eral relativistic spin-spin interactions in Hawking radiation from rotating black holes, 



see 



De-Chang and Stoikovid (120101 ). 



Great efforts have been devoted so far to work out, with a variety of techniques, several 
kinds of next-to-leading order spin-spin and spin-orbit effects in compact binary systems 
for various ranges of masses, spins and orbital geometries in view of the impac t that they 
would have on the emission of gravitational waves; see the review by iBlanchetl (j2006[ ). We 
recall that, actually, several gravitational waves observatories are planned- and also already 



opera ting-to view binaries in various stages of their evolution ( lGiazottdl2008l : iFairhurst et al. 
201ll ). so tha t they might detect them. For two-b o dy systems eff e cts in v olving sp ins 



linearl^, se^ Konigsdorffer and Gopakumar ( 2005 ): D amour et al. ( 2008 ): Tessmeij ( 2009 ): 
Tessrner et al.l (120101): ICornish and Keyl (120101 ) . See also ICergelv et al.l (ll998al JblJcr): iGergely 
(1l999l . l2010h : iMaiar and VasuthI (l2006l . l2008h : IVasuth and Maiai] (120071 ) in which, among 
other things, d i fferen t sets of independ e nt va riables were used. We also mention that 
Steinhoff et al.l ( 120081 ) and iHergt et al.l ( 120101 ) investi gated next-to-leading-order s pin- 
squared dynamics of general compact binaries, while iHartung and Steinhofla ( 120111 ) dealt 
with the next-to-leading order spin-orbit and spin-spin effects for a system of N gravitating 
spinning compact objects. In general, we will not deal in more details with such a host of 
gravitational wave-related features of motion since, as it will be clear a posteriori from our 
analysis itself, the current observational methods are quite insufficient to resolve such very 
peculiar spin-spin and spin-orbit effects. A direct and explicit comparison of our results 
with some of the leading order ones existing in literature would be difficult to be made in 
a meaningful way because of the differences in the approaches and the languages followed, 
driven by different original tasks. The tone of our work is mainly phenomenological, so that 
we choose the most widely adopted parameterization in interpreting the empirical analyses 
usually performed in astronomical and astrophysical scenarios like the planetary motions in 
the solar system and binary pulsars. 



^We mean that each spin enters the equations with the first power. 



^For example, the results by iKonigsdorffer and Gopakumar! (l2005f ) are valid for binary 
systems w ith equal masses and arb itrary spins, and arbitrary masses with only one of them 
spinning. 



Cornish and Keyl (120101 ) derive a complete analytic solution for binary motion 



with spin-orbit interaction. 
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The plan of the paper is as follows. In Section |2] we analytically work out the long-term 
temporal changes of order (9(c~^) for all the Keplerian orbital elements of spinning and 
non-spinning particles orbiting a massive central source. Contrary to some other works 
on such topic existing in literature, we will not restrict ourselves to specific, simplified 
configurations: our results will be quite general, as far as both the spatial orientations of 
the spins and the orbital geometry of the particle are concerned. In Section [3] we apply our 
results to the aforementioned specific astronomical and astrophysical scenarios. Section H] 
summarizes our findings. 



2. Analytical calculation 

In order to evaluate the long-term effects caused by Eq. (fl]) on the orbital motion 
of the moving particle, it must, first, be evaluated onto the unperturbeci^ Keplerian 
orbit and averaged over one orbital revolution. A similar approach can be found in 



Damour and Deruelld ( 1l985l ). The orbital motion in the two-body problem can be usefully 



param eterized in terms of the standard osculating Keplerian orbital parameters (jCapderou 



20051 ). They are the semimajor axis a and the eccentricity e, which fix the size and the 
shap^ of the ellipse, and the longitude of the ascending node i?, the argument of pericenter 
u and the inclination J, which determine the orientation of the ellipse in the inertial spac^. 
The instantaneous position of the test particle along the Keplerian ellipse is reckoned by 
the true anomaly /, which, in combination with u, yields the argument of the latitude 
u = u + f. The unperturbed Keplerian mean motion is n = a^; it is connected 

with the mean anomaly Ai through Ai = n{t — tp), where tp is the time of passage at 
the pericenter. The Keplerian mean motion is also related to the orbital period Pb by 
n = 27r/Pb- 



^Notice that, in princ i ple, it would be, perhaps, possible to use the perturbative techniques 
devised by ICalura et al.l ( 119971 . Il998l ) for various first order post-Newtonian reference orbits 
involving non-spinning bodies, but they would yield additional negligible mixed orbital effects 
of order higher than (9(c~^). 

^It can be < e < 1: for e = the ellipse reduces to a circle. 

''They can be thought as the three Euler angles establishing the fixed orientation of a 
rigid body, i.e. the Keplerian ellipse which changes neither its size nor its shape, with respect 
to an inertial reference frame. 
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By recalling that on the unperturbed Keplerian ellipse it is fICapderoul |2005| ) 



X 



< y 



dt 



a(l-e^) 
1+e cos / ' 

r (cos J? COS u — COS / sin i? sin u) 
r (sin i? cos u + cos / cos J? sin u) 
r sin /sin u, 

n(l+ecos/)2 ' 



(2) 



it turns out 



2c2a3(l-e2)3/2 



{2 (cr • J) + 3 ( Ja; cos f} + Jy sin i?) ((Tj^ cos i? + cr^ sin i?) + 



+ [Jz sin / + cos / {Jy cos Q — Jx sin i?)] [a^ sin / + cos / {cjy cos Q — (Jx sin i?)]} , 

(3) 

in which the brackets (■ ■ ■ ) denote the average of their content over one orbital period. We 
computed it by using the true anomaly as independent variable: 



n 



(4) 



'0 ^" Jo 

by means of Eq. Notice that, in principle, other averaging procedures implying the 
use of the mean anomaly Ai and the eccentric ano r naly E as independent variables may 
have been adopted; see, e.g., iGopakumar and lyerl (119971 ) for one based on the use of 
E. It turns out that, in the present case, the use of / through Eq. ([2]) allows to obtain 
exact expressions in e. Instead, it would not be possible with Ai and E since necessarily 
approximatecl§ expressions for cos / and sin / or, equivalently, of cos E a nd sin-E, wh ich 



enter 'Ha,'HnT,,'H .iT,, in terms of M. to some order in e should be used (jSmartI Il961 



KovalevskyI Il963l ) . The dimensions of {'Ha) are correctly [{'Ha)] = T~^. In the case 
of 'HaL, it must be recalled that, in the unperturbed case, the reduced orbital angular 
momentu m L is a ve ctor directed along the out-of-plane direction, determined by the unit 
vector n (jSoffellll989h 

sin/ sin i?, 



n = < — sin /cos i?, 



(5) 



cos /, 



'^It is, e.g., cos / = cos - e + ecos27\/i + C(e^), sin/ = sin + esin2A^ + 0{e'^). 
For cosE and sinE as inf inite summations of trigonometric functions over Ai, see, e.g.. 



Tessmer and Schaferl (120111 ). 
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and whose magnitude is (ICapderoull2005l ) 

L = 

Thus, 



na' 



SGMn 
2c^a(l - e^' 



[a^ cos / + sin I (a^ sin Q ~ ay cos i?)] : 



(6) 



(7) 



also in this case, [("Ho-l)] = T~^, as expected. Let us note that the results of Eq. ([3]) and 
Eq. (jTj) are exact in the sense that approximations neither in e nor in / have been used to 
obtain them. Moreover, Eq. ([3]) and Eq. ([7]) are defined also for e — t- 0, / — 0. 

It is, now, possible to straightforwardly work out the secular variations of all the 
six Keplerian orbita l elements through the Lagrange planetary equations. They ar^ 
fiBertotti et al.ll2003h 



/■ da 

dt 

de 
dt 



dl 

dt 



j_ an \ 

' na \ dM) 



1-e" 



1 dn 



dM J ' 



I \ an duj ) 



dn 

dt 



dm 
dt 



dM 
dt 



1 rdTl\ 

^^sinl \ dl ) 



an I tan(j/2) an 
e ae ^ ai 



n + 



2a^ 

aa 



an 

de 



in which w = i? + is the longitude of the pericentre, one of the parameters usually 
estimated by the astronomers when they process long data records of the planets of our 
solar system, and TZ denotes the averaged perturbation of the Newtonian gravitational 
potential: in our case it is 

7^ = (H.) + {H,l) ■ (9) 

A rapid inspection of Eq. (|3]) and Eq. ([7]) tells us that, according to Eq. ([H]), all the 
Keplerian orbital elements should, in principle, experience secular variations, apart from 
the semi-major axis a. 



^Note that the convention Ug = GM/r, yielding F = ^Ug, is often used in many 
textbooks for the Newtonian gravitational potential Ug] in such a case, the right-hand-sides 
of the Lagrange planetary equations have a minus sign with respect to Eq. ([8]). 



dt 



2.1. The spin-spin Stern-Gerlach precessions 

Actually, for of Eq. the secular changes are 

= 0, 
— = 

W = 2^3^5|^3^ {sin J [(Jy(Ta; + JxCTy) COS 2i? + (J^a^ - J^a^;) sin 2i7] - 

— cos / [{JzCTx + JxO'z) COS i? + ( JzCTy + JyO"^) siu i?] } , 
^ = 2e2a5n(l-e2)2 {^OS 2/ CSC / [{Jz<Jx + ^xCT^) siu i? - ( J^CTj^ + Jj^CT^) COS J7] + 
+ COS I [J • Cr — SJzCTz + ( JyCTy — JxCx) COS 2 J? — ( JxCy + JyO'x) siu 2 J?]} , 

^ = 2^2^5|^3^ {2J • cr - 3 (J^. COS i? + sin i?) (cr^ COS i? + CTy sin i?) - 



— 3 [Jz sin / + cos / {Jy cos Q — Jx sin J?)] [a^ sin / + cos / {ay cos Q — ax sin i?)] — 

— tan {1/2) [Jz cos / + sin / {—Jy cos i? + sin i?)] [a^ sin J + cos / {ay cos Q ~ ax sin i?)] - 

— tan {1/2) [Jz sin / + cos / {Jy cos Q — J^ sin i?)] [cr^ cos / + sin I {—ay cos Q + a^ sin i?)]} 

^ = ~ 2c2a5na^e2)3/2 ( J • O") + 3 (J^; COS Q + Jy siu i?) (cTa; COS i? + (7^ siu i?) + 

+ 3 [Jz sin / + COS / {Jy COS f2 — J^ sin i?)] [a^ sin I + cos / cos Q — a^ sin i?)] } . 

(10) 

Their dimensions are correctly those of T ^. 
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2.2. The spin-orbit precessions 

2.2.1. The case of a small spinning object orbiting around a non-rotating central body 
The precessions caused by "Hq-l in Eq. ([7]) are simpler, amounting to 
= 0, 
= 0, 



da 
dt 



de 
dt 



dl 3GM{ax cos H+Uy sin J7) 

dt ~ 2c2a3(l-e2)3/2 ' 

(11) 

dn 3GA/[(72 +cot /(cTy cos J7—(T3; sin J7)] 

'dt' ~ 2c2a3{l-e2)3/2 ' 

d-uo 3GM{2[az cos I+s'mI{ax sin fl—cTy cos I?)] — [ctz sin7+cos I{c!y cos fl—a^ sin J7)] tan(//2)} 

"dT ~ 2c2a3(l-e2)3/2 ) 

AM ^ n. 

also in this case, it is easily recognized that their dimensions are those of T~^, as expected. 
The results of Eq. (fTOj) and Eq. (ITT]) are completely general, and are exact both in e and 
in /. Moreover, they are valid also for e ^ 0, / — )■ 0, apart from the node precessions 
which become singu l ar for / = 0. They may be compared with those obtained by 
Barker fc O'Connell fll97nh: concerning the occurrence of the change in the inclination, it 



was predicted by 



Barker fc O'Connelll (flOToh . 



2.2.2. The case of a non- spinning particle orbiting around a rotating central body: a 

generalization of the Lense-Thirring effect 



Another larger term should, actually, be considere d in addition to 'Hrrr. in E g- dl])- It 
comes from the full two-body spin-orbit Hamiltonian ( iBarker &: O'Connelll 119751 ): in the 
extreme mass ratio limit M 3> m, it can straightforwardly be obtained from Eq. ([1]) with 
the replacement 

-Ma 2J. (12) 



Thus, we have 



n 



JL 



2G 

^2 



iJ-L) 



(13) 
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it is \HjL__^ 

^ jloriolboOll ) 



terr 



T ^. It can easily be recognized as just the gravitomagnetic Lense-Thirring 



JL 



V 



gm 



with the gravitomagnetic vector potential given by (IMashhoonl 120071 ) 

_ 2GJ X r 



(14) 



(15) 



[Ag] = T~^. Then, Eq. (fTTI) . with the replacement of Eq. (fT2|) . straightforwardly gives 
us the generalized Lense-Thirring precessions of the Keplerian orbital elements of a test 
particle for an arbitrary orientation of the proper angular momentum of the source. They 
are 



f da 

dt 

de 
dt 

dJ 
dt 

dn 

dt 



dzu 
dt 

dM 



0, 



0, 



2G(Jx cos Q+Jy sin Q) 
c2a3(l-e2)3/2 ; 

2G[Jz+cot I{Jy COS H—Jx sin H)] 
c2a3(l-e2)3/2 ' 

2G{2[Jz cos /+sin/(Jx sin Q—Jy cos Q)] — {Jz sin/+cos I{Jy COS ^? — Jx 

sin J?)] tan(//2)} 
c2a3(l-e2)3/2 ' 



(16) 



dAa n. 

thei r dimensions are co r rectly those of T~^. On the contrary, all the usual derivations of 
the iLense and ThirringI (119181 ) precessions of the Keplerian orbital elements existing in 
literature are based on the particular choice of aligning J along the z axis. To this aim, 
it can be noted that Eq. yieldfl] just the usual Lense-Thirring rates for J^, Jy — )■ 0. 
Eq. f ll6l) tells us that also the inclination I experiences a secular gravitomagnetic change 
if J has components along the x and y axes as well: it is independent of the inclination / 
itself. Moreover, and Jy induce additional precessions for the node i? and the longitude 
of pericentre w which depend on /. In general, all the additional precessions depend on J? 
as well. 



^''See also iLandau and Lifshitzl (119751 ) and lAshby fc AllisonI (119931 ). 

^^Indeed, by posing 1/2 = a, it is easy to show that — 2cos/-|-sin/tan(//2) = 1 — 3 cos/. 
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3. Observability of the computed effects 

In this Section, we compare our predicted results to the latest observationally 
determined orbital precessions in some astronomical and astrophysical scenarios. We will 
work them in detail, also because the results obtained here for the spin configurations may 
turn out to be useful also for different dynamical effects. 



3.1. The Sun-Mercury system 

Let us, first, apply our results to the solar system. In this case, the inclination / may 
typically be referred, e.g., to the ecliptic plane, assumed as reference {xy} plane, and the 
longitude of the ascending node i? is counted from the Vernal equinox, assumed directed 
along the positive x axis. Concerning the Sun's spin axis, its orientation with respect to the 
J2000 ecliptic is 



T0 



J? 



cos 5^ cosa^ = 0.122, 

cos(5|f sinoQ cose + sin(5[f sine = —0.031, 

— cos 5q sin sin e + sin 5^ cos e = 0.992, 



(17) 



where 



an 



286.13 deg. 



63.87 deg, 

are thj^ right ascension and declination of the Sun's north pole of rotation, respectively 
( ISeidelmann et al.l 120071 ). while 

e = 23.439 deg (19) 



is the obliquity of the Earth's equator to the ecliptic at J2000 fiFukushimal l2003l ) . The 
magnitude of the Sun's pr oper angular momentum was measured with the helioseismology 
technique. It amounts to (jPijpersI 119981 ) 



Jq = (190.0 ± 1.5) X 10^^ kg m^ s'^ 



(20) 



^^The right ascension a and the declination S refer to the mean terrestrial equator at 
J2000. 
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Instead, for Mercury we have 



A? 5 
COS Oq cos 



0.091, 



a. 



from ( ISeidelmann et al.l 120071 ) 



: COS 5^ sin cose + sin ^ sine = —0.081, 
cos (5^ sin sin e + sin 5^ cos e = 0.993, 



al = 281.01 deg, 
= 61.45 deg. 



(21) 



(22) 



and 



a 



5 



2.4 X 10^ s"\ 



(23) 



as can be inferred from the values of its equatorial radius, angular rotation spee d and 
normalized polar moment of inertia!^ ( Tholen et al. 2000 : Seidelmann et al. 2007 ). 



Since ( iMurray and Dermottlll999l ) 



0.387 AU, 
0.205, 
7.005 deg, 
48.330 deg. 



(24) 



it turns out that the spin-spin and spin-orbit precessions of Mercury are up to 5 — 11 orde rs 
of ma gnitude smaller than the present-day level of accuracy in measuring them fjPitjeva 
2010l ). i.e. about a few milliarcseconds per century (mas cty~^ in the following). Indeed, 
Eq. ( JTOj) yields for the "Ho-— related precessions 

^ -4 X 10-1^ mas cty-\ 



dt 



dt 



—5 X 10 mas cty ^, 
: 1.35 X 10"^ mas cty~\ 
dM. = 1.33 X 10-9 mas cty-\ 



dm 
dt 



(25) 



13 



See also |http : / / nssdc . gsf c . nasa. gov / planetary / f actsheet / mercury fact . html on the WEB 
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while the Mercury-induced spin-orbit rates of Eq. f[TT|) are 

: 3 X 10"^ mas cty~^, 



( dl 

dt 



dt 



dvj 
dt 



-1.2 X 10"^ mas cty"\ 



-3.8587 X 10"^ mas cty" 



(26) 



From the extremely small values of the secular rates of the inclination it can be inferred 
that the general relativistic spin-spin and spin-orbit effects considered here did not play any 
role in the solar system's evolution throughout its lifetime of about 4.5 Gyr. 

Moving to the Lense-Thirring precessions of Eq. (fT6l) . they are 

: 0.06 mas cty^^, 



dl 

dt 



dQ 
dt 

dzj 
dt 



0.08 mas cty" 



-2.01 mas cty 



(27) 



The secular rate of the inclination of the orbit of Mercury is, at present, too small by likely 
two orders of magnitude to be measurable in the near future. An important result for the 
node is that the additional precession due to J® and J®, as large as —0.92 mas cty~^ 



I — I r 

almo st cancels the precession caused only by Jz, which amounts tcLJ 1 mas cty~^ (jlorio 



20051 ). Thus, the net Lense-Thirring node precession is two orders of magnitude smaller 
than the present-day accuracy in determining the nodal rates for the inner planets. On 
the contrary, the total precession of is pract ically equal to the valuJ^ —1.98 mas cty"^ 
obtained by assuming = Jy = (llorid 120071 ). A recent application of such results to the 
motion of Mercury in vie w of the late st observations from the MESSENGER spacecraft 
orbiting it can be found in llorid (120111 ) . 



Finall y, we point out that the still unmodelled (and undetected) 2PN pointlike perihelion 



precession (jPamour and Schaferlll988l : IWexlll995l ). proportional to 'dj2PN (GM)^'^c 



-7/2 



up to terms of order O(e^), is about 10 ^ mas cty ^ for Mercury. 



3.2. The double pulsar PSR J0737-3039A/B system 



A relatively more promising sc e nario is represen ted by the double pulsar PSR 
J0737-3039A/B (IBurgay et al.l l2003l : iLyne et al.ll2004l ): it is a close binary system at 



^ iCugusi and Proverbid (119781 ) gave a value too large by one order of magnitude. 



15 



Cugusi and Proverbid (11978) andlSoffell (119891) re l eased values larger by a factor 10 — 5 



respectively, while the figure by Barker fc O'Connelll (jl970[ ) was only 1.5 times larger. 
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aboutcj 500 pc from us ( Kramer et al. 2006 ) constituted of two neutron stars both visible 
as radioDulsara^j orbit ing along moderately eccentric barycentric ellipses in 2.45 hr= 0.102 



d (IKramer et al.ll2006l ). In this case the inclination / refers to the plane of the sky, assumed 
as reference {xy} plane, so that the z axis is directed towards the ob server along the 
line-of-sight direction. The relevant parameters of such a system arj^ ( Kramer et al. 2006) 



' m 


= Ma + Mb = 2.58708 M 


Ma 


= 1.3381 Mq, 


Mb 


= 1.2489 Mq, 


r 

u 


= ^ = 1.0714, 

Mb ' 

= „fAA/B ^ 0.6459 M„, 
= 44 05 Hz 


< I^B 


= 0.36 Hz, 


a 


= 8.78949386 x 10^ m. 


e 


= 0.0877775, 


I 


= 88.69 deg. 


fiL 


= 7.03066 X 10^^ kg m^ s" 


Ja 


= 2.768 X 10^° kg m^ s'^, 




= 2.3 X 10^8 kg m2 s-\ 



(28) 



^^Actually, doubts on the reliability of such an estimate may be cast since it is based on a 
model of the interstellar electron den sity. On the other hand, subsequent very long baseline 
interferometry (VLBI) observations (IDeller et al.ll2009l ) yield a distance of 1150l^go P*^- 

^''Actually, B disappeared from our view in 2009 (jPerera et al.ll2010l ). 

^'^Concerning the inclination J, in the PSR J0737-3039A/B system it is, actually, inferred 
from the determined value of sin J through the Shapiro delay measurements. Thus, in ad- 
dition to the quoted value of J = 88.69 deg, also the value I = 180 deg— 88.69 deg= 91.31 
deg is, in principle, admissible. From a practical, calculational point of view, this has no 
appreciable consequences on the numerical results presented below. 
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where we used the standard value of the pulsar's moment of inertia X !^ 10^^ kg m^ and the 
fact that the rotational s pin periods of A a nd B are Ta = 1/z^a = 23 ms (IBurgay et al.ll2003l ) 
and Tb = I/i^b = 2.8 s (ILyne et al.ll2004l ). respectively. From Eq. ( l28l) it can be noticed 



that the longitude of the ascending node i? does not figure in the list of known parameters 
of PSR J0737-3039A/B. It is so because Q does not directly affect the direct observable 
quan tity, i.e. the standard timing formula (jPamour and Deruelldll986l : iTaylor and Weisberg 
19891 ). so that it cannot, i n gener a l, be determined from pulsars' timing data analysis. For 
a general review, see, e.g., IStaird (120031 ). The node could be measured only if a pulsar is 
close enough to the Earth. Indeed, in this case the orbital motion of the Earth changes 
the apparent inclination angle / of the pulsar orbit on the sky, an effect known as the 
annual-orbital parallax flKopeikiij Il995l ) . It results in a periodic change of the projected 



semi-ma 



lor axis. A seco nd contribution comes from the transverse motion in the plane of 
the sky fjKopeikinlll996l ). yielding a secular variation of the projected semi-major axis. By 



including both these effects in the model of the pulse arrival times, the longitude of the 



ascen ding node i? can be determined, as in the case of PSR J0437-4715 (Ivan Straten et al. 
200ll ). which is at only 140 pc from us. 



Concerning the orientations of the spin axes, it turns o ut that J a is substantia lly 
aligned^ with the t otal angular momentum of the system ( Manchester et al. 2005 : 
Ferdman et al.ll2008l ). which practically coincides with the orbital angular momentum /xX, 
as can be inferred from Eq. (I28l) . Thus, according to Eq. ([5]), we can pose 



J. 



J 



0.997 sin J7, 



-0.997 cos J7, 



(29) 



0.023 



Ja lies almost entirely in the plane of the sky. Instead, ctp which undergoes th e general 
relativistic de Sitter precessioio ( iBarker fc O'Connell 1975 : Boerner et al. 1975 ) descr ibing 
a full cycle in 75 yr, has a different orientation (IBreton et al.l 120081 : iPerera et al.ll2010[ ). At 



^^Actually, according to iFerdman et al.l ( 120081 ). there is still a 30% chance that the tilt of 
A's spin is larger than 6 deg. This implies that the projection of Ja on the orbital plane 
may be larger than 10% of the total. 

^°It is the main candidate to explain the disappearance of the radiopulses of B in 2009; if 
it is the sole cause of such a phenomenon, t he B's signal should reappear in about 2035 or 
even earlier, depending on the beam shape (jKramerll2010l ). 
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the epoch May 2, 2006 (MJD 53857), according to Figure 1 of Breton et all f l2008f ). it is 

^ sin (/) cos 6^ = -0.501, 



< 



a 



y 



where 



sin sin 6' = 0.597, 
cos = 0.625, 

= 51.21 deg, 
= 130.02 deg. 



(30) 



(31) 



Actually, the frame used by iBreton et al.l (|2008[ ). dubbed {x\y\z'} by us, is different from 
the one {x, y, z} used here. Indeed, the x axis is directed along the line-of-sight towards 
the Earth, the z axis is directed along the (projected) orbital angular momentum, i.e. it is 
z = sin In, and the y axis coincides with the line of the nodes. In other words, the y and 
the z axes lie in the plane of the sky, being rotated by i? with respect to the x and y axes. 
Thus, a^, are linear combinationcj of a^, and a^, involving i? as well, while = a^,. 

Concerning the spin-spin and spin-orbit Hamiltonians of Eq. ([T]), they are to be 
modified in order to take into account the fact tha t the mass of B is almost equal to 
that of A. According to f lBarker fc 0'Connelllll975[ ). Ha must simply be re-scaled by the 
dimensionless factor 

Mb 



1.933, 



(32) 



while in the spin-orbit term of Eq. ([T]) the replacement 



1 + 



A ^ B 



3 Ma 

must be done. Moreover, a second term containing the spin of A must be added, so that 

3 



(33) 



G 




4Mb\ 




_^Ma(h 


^2(^ 3 




3 Ma; 



-Mf 



(34) 



here jt a = Ja/Ma. In order to facilit ate a comparison with other works existing in literature 
like, e.g., ID amour and Taylorl (jl992[ ) it is convenient to re- write Eq. fjM|) . after a little 
algebra, as 



JL 



2G 

^2 ^ 3 



1 + 



3 Ma 
4Mr 



L)+ 1 + 



3 Mr 



4Ma 



(Ja-L) 



(35) 



^^They are = cos i?(T^, + sin^ / sin Ha^, and = sin i?(T^, — sin^ I cos Ha^, . 

■J 11 z y V z 
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Thus, the resulting orbital effects can straightforwardly be inferred from Eq. (1161) with the 
replacement 

(36) 



1 



4 Ma 



for the precession by Ja and vice-v ersa for the effects induced by Sb- Some of them can be 
compared with those worked out by lDamour and Taylon (jl992l ). For example, our expression 
for the precession of / due to, say, Ja obtainable from Eg. (fT6l) with the substituti on of Eq. 
f p6|) coincides with Eq. (3.27), written for Eq. (3.7), of lDamour and Tavlorl (119921). Indee d. 
the un it vector I = i of the line of the nodes entering Eq. (3.27) of lDamour and Taylor 
(1992) has cartesian component^ 1^ = 1^. = c os i?, I„ = i„ = sin i? , so th at its scalar 
product with the the spin of A in Eq. (3.27) by ID amour and Taylorl (Il992[ ) yields just the 
geometrical factor of our Eq. ( ITBl) . 

As far as the spin-obit precessions caused by Ja are concerned, Eq. ([5]) and Eq. ( IT6|l 
tell us that, since we are assuming the spin axis of A exactly aligned with the orbital 
angular momentum, the precessions of / and Q vanish, contrary to that of u which becomes 
proportional to cos^ / + sin^ /: the dependence on Q disappears. A similar situation occurs 
for the Sun-Mercury system in which, as we have seen in Section 13. H the Lense-Thirring 
precessions of / and Q are two orders of magnitude smaller than the perihelion precession. 
The magnitude of the periastron precession due to Ja is of the order of 



duj 

r 

dt ' 



3.7 X 10"^ deg yr"^; 



(37) 



the present-day level of accuracy in determining t he periastron preces sion of the double 
pulsar from observations is 6.8 x 10""^ deg yr~^ ( Kramer et al.l 120061 ). The issue of 



investigating t he actual rn e asurability of the spin-or bit periastron precession caused by Ja 
was tackled in llorid (120091 ): iKramer and Wexl (|20C)9|). along with the aliasing effects of t he 
precessional effects caused by the IPN and 2PN (jPamour and Schaferlll988l : IWexlll995l ) 
non-spinning point particles contributions. Here we briefly recall that the nominal value of 
the 2PN precession, not yet measured, is of the order of 10~^ deg yr~^, with an uncertainty 
of 10 ~^ deg yr~ ^ due to our imperfect knowledge of the system's total mass and semi-major 



axis (lloridl2009l ). Instead, the uncertainty in the measured IPN prece ssion due to the errors 
in the system's parameters is of about 3 x 10~^ deg yr~^ (loriol 20091). If one tries to take 



into account a possible tilt of Ja with respect to n (IFerdman et al.ll2008l ). the effects would 
be very small and undetectable. Indeed, Eq. (fT6|l . with Eq. (136|1 . tells us that the non-zero 
secular variation of, say, the inclination would be 



dl 



= jA cos + sin 1.8 X 10"^ deg yr 



(38) 



22c 



See Eq. (2.12a) of lDamour and Tavloi] (Il992l ). 
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where it is intended that the components of Ja do not coincide with those of n given in 
Eq. (j5]). Furthermore, dl/dt is not an observable quantity in the timing observations of 
the double pulsar. Indeed, they are sensitive to a change in the projected semi- major axis 
( jPamour and Taylorlll992l ). which is proportional to sin J. This fact, on the other hand, 
supports the previous argument that the effect on I is negligible from an observational point 
of view since the rate of variation of the projected semi-major axis is proportional to cos/, 
which is small for the pulsar PSR J0737-3039A/B. Instead, all the spin-orbit precessions 
due to Sb do not, in general, vanish: they are as large as up to 1 — 2 x 10"^ deg yr~^. The 
spin-spin Stern-Gerlach rates are quite smaller, being all of the order of 2 — 4 x 10"^^ deg 
yr~^ 



4. Summary and conclusions 

In this paper, we analytically worked out the secular variations of the six osculating 
Keplerian orbital elements of a spinning particle in free motion around a central, slowly 
rotating body caused by the general relativistic spin-spin and spin-orbit interactions of 
order 0{c~'^) in the weak-field and slow-motion approximation. We neither restricted to 
specific spatial orientations of the angular momenta involved nor to particular orbital 
configurations: the results obtained are, thus, exact with respect to such issues and quite 
general. We adopted the Lagrange planetary equations after having averaged the perturbing 
reduced Hamiltionians over one orbital revolution. As reference, unperturbed orbit we 
adopted the Keplerian ellipse: far smaller effects due to mixed static-spin general relativistic 
terms could have, in principle, been worked out, but we neglected them because of their 
negligible magnitude. 

Concerning the spin-spin Stern-Gerlach-like term, it turns out that the inclination, 
the node, the longitude of pericenter and the mean anomaly undergo secular precessions 
depending on the particle's semi-major axis as a~^n~^ oc a~^/^ and the eccentricity. The 
dependence on the particle's inclination and node is rather complicated; all the components 
of both spins enter the formulas as well. 

The spin-orbit term containing the particle's spin causes secular precessions of the 
inclination, the node and the longitude of pericenter which go as a"^. The spin-orbit term 
containing the angular momentum of the central body, arising from the full two-body 
spin-orbit Hamiltonian, yields the usual Lense-Thirring effect for a specific orientation of the 
body's spin. Instead, for a generic orientation of it we obtained secular precessions for the 
inclination, the node and the longitude of pericenter which depend on all the components 
of the angular momentum of the central body, and on the node and the inclination of the 
orbital plane. The dependence on the semi-major axis and the eccentricity is as for the 
other spin-orbit term. 

We generalized the results obtained to a generic two-body system with arbitrary 
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masses. 

We applied our predictions to some specific astronomical scenarios like the Sun- Mercury 
system and the double pulsar. Concerning Mercury, it turns out that the Stern-Gerlach 
spin-spin precessions are too small to be detected by 9 — 1 1 orders of magnitude with respect 
to the present-day level of accuracy in empirically determining the planetary precessions, 
which is a few milliarcseconds per century. The precessions caused by the Mercury's spin are 
too small by 5 — 11 orders of magnitude. In regard to the Lense-Thirring effect, the slight 
misalignment of the Sun's equator with respect to mean ecliptic at J2000, usually adopted 
in all the planetary data reductions performed by the most influential teams of astronomers, 
induces a tiny secular precession of the inclination of about 0.06 milliarcseconds per century. 
Concerning the node, its secular precession turns out to be reduced to 0.08 milliarcseconds 
per year with respect to the typical value of 1 milliarcseconds per year usually computed 
by neglecting the departure of the solar equator from the ecliptic, while the precession of 
the longitude of the perihelion is slightly augmented to —2.01 milliarcseconds per year with 
respect to the standard value of —1.98 milliarcseconds per year. Such results are important 
in view of a possible reliable detection of the solar Lense-Thirring effect in the near future. 

Moving to the double pulsar scenario, the Stern-Gerlach spin-spin precessions are 
about 8 orders of magnitude smaller than the present-day accuracy of 6.4 x 10~^ degrees 
per year in determining the orbital precessions from timing data. The spin-orbit effects due 
to the angular momentum of B, which has the slowest rotation in the pair, are too small 
by 2 orders of magnitude. Instead, the spin-orbit precessions caused by the spin of A are 
all vanishing because of its alignment with the system's orbital angular momentum, apart 
from the periastron which precesses at a rate of —3.7 x 10~^ degrees per year. 

The full generality of our calculations with respect to the orbital and the spin spatial 
configurations allows them to be straightforwardly extended to other systems like extrasolar 
planets. Indeed, they present a rich variety of orbital configurations and spins orientations 
which make them potentially interesting candidates to apply our results, especially those 
regarding the generalized Lense-Thirring precessions. This may be the subject of a 
forthcoming paper. 
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